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The Intrinsic Alignment of Dark Halo Substructures 

JOUNGHUN LEE^’^, XI KANG^, AND YiPENG JlNG^ 

ABSTRACT 

We investigate the intrinsic alignments of dark halo snbstrnctnres with their 
host halo major-axis orientations both analytically and nnmerically. Analytically, 
we derive the probability density distribution of the angles between the minor 
axes of the substructures and the major axes of their host halos from the physi¬ 
cal principles, under the assumption that the substructure alignment on galaxy 
scale is a consequence of the tidal helds of the host halo gravitational potential. 
Numerically, we use a sample of four cluster-scale halos and their galaxy-scale 
substructures from recent high-resolution N-body simulations to measure the 
probability density distribution. We compare the numerical distribution with 
the analytic prediction, and find that the two results agree with each other very 
well. We conclude that our analytic model provides a quantitative physical ex¬ 
planation for the intrinsic alignment of dark halo substructures. We also discuss 
the possibility of discriminating our model from the anisotropic infall scenario by 
testing it against very large N-body simulations in the future. 

Subject headings: cosmology:theory — large-scale structure of universe 


1. INTRODUCTION 

The dark halo substructure has recently come to one of the most lively topics in cosmol¬ 
ogy. Although the standard cosmological paradigm based on the cold dark matter (CDM) 
concept generically predicts the presence of the substructure inside the dark matter halos, 
there are many questions yet to be answered associated with the dark halo substructures. 
The intrinsic alignment effect of the dark halo substructure is one of those questions. 
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There are plenty of observational evidences that the major axes of the brightest cluster 
galaxies (BCGs) have a strong tendency to be aligned with that of their host clusters (Sastry 
1968; Carter & Metcalfe 1980; Binggeli 1982; Struble & Peebles 1985; Rhee & Katgert 1987; 
West 1989, 1994; Plionis 1994; Fuller, West & Bridges 1999; Kim et ah 2002). The most 
popular theory for the BCG alignment is the anisotropic infall scenario based on the standard 
hierarchical clustering model (West 1989): The initial density held of CDM is web-like, 
interconnected by the primordial hlaments (Bond 1987; Bond, Kofman, & Pogosyan 1996). 
The gravitational collapse and merging to form structures occurs not in an isotropic way but 
in an anisotropic way along the large-scale hlaments. Accordingly, the infall of materials into 
a cluster also occurs along the primordial hlament, which will induce the alignment between 
the orientation of a host cluster and that of BCG embedded in it. 

There are several reasons that the anisotropic infall theory became so popular: Being 
simple and intuitive, it hts very well into the cold dark matter paradigm. In addition, 
it has been supported by several numerical simulations (e.g.. West, Villumsen, & Dekel 
1991; van Haarlem & van de Weygaert 1993; Dubinski 1998; Faltenbacher et al. 2002) which 
demonstrated that the gravitational infall and merging of materials indeed occurs along the 
hlaments. 

Nevertheless, the theory is only qualitative and still incomplete. Recent observations 
indicate that not only the BCGs but also the less dominant cluster galaxies exhibit the 
alignment ehect to a non-negligible degree (Plionis & Basilakos 2002; Plionis et al. 2003; 
Pereira & Kuhn 2004). In the anisotropic infall model, the substructure alignment is a 
primordial ehect, and would get damped away quickly by the subsequent nonlinear processes 
such as the violent relaxation, the secondary infall, and so on (Quinn & Binney 1992; Coutts 
1996). Therefore, it is very unlikely that the cluster galaxies other than the BCGs keep the 
primordial alignment ehect till the present epoch (Plionis et al. 2003). Here, we propose 
that the initial tidal interaction between the subhalos and the host halo will be responsible 
for the observed intrinsic alignment of the cluster galaxies. 


2. ANALYTICAL PREDICTIONS 


When a subhalo forms inside a host halo, it acquires the angular momentum L = (LQ 
due to the tidal shear held T = (T^) generated by the gravitational potential of the host 
halo (T): Tij = didj^. Lee & Pen (2000, 2001) proposed the following formula to quantify 
the mutual correlations between T and L: 


{LjL,|T> = 


1 + c 


5ij cTif^Tf^j. 
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where c G [0,1] is a correlation parameter to quantify the strength of the correlation between 
T and L, and T = (Tjj) is a unit traceless tidal shear tensor dehned as Tij = Tjj/|T| with 
Tij = Tij — Tr(T)5jj/3, and L = (Lj) is a rescaled but not a unit angular momentum. If 
we replace the rescaled angular momentum by the unit angular momentum, Lj = L/|L| in 
equation (1), then the correlation parameter c is reduced by a factor of 3/5. Note here that 
the LHS of equation (1) represents a conditional ensemble average of LiLj provided that the 
unit traceless tidal shear tensor is given as Tij. For the detailed explanations of equation 
(1), see Appendix A in Lee & Pen (2001). 

It is naturally expected that c decreases with time as the correlation between T and L 
must decrease after the moment of the turn-around due to the subsequent nonlinear process. 
Lee & Pen (2002) found c ~ 0.3 at present epoch by analyzing the data from the Tully 
Galaxy Catalog and the Point Source Redshift Catalog Redshift Survey (in their original 
work, they used a reduced correlation parameter a = 3c/5 and found a = 0.18). 

Strictly speaking, equation (1) holds only if T and L are defined at the same positions 
(Lee & Pen 2000, 2001). Here, they don’t; T and L are dehned at the centers of the mass of 
the host halo and the subhalo, respectively. For simplicity, here we just assume that equation 
(1) still holds, ignoring the separation between the centers of the mass of the subhalos and 
that of the host halo. 


The distribution of L under the inhuence of the tidal held is often regarded as Gaussian 
(Catelan & Theuns 1996; Lee & Pen 2001): 


P(L) 


1 

[(27r)3det(M)]V2 




( 2 ) 


where the covariance matrix Mij = (LjLjjT) is related to the tidal shear held T by equation 
(1). In the principal axis frame of T, let us express L in terms of the spherical coordinates: 
L = (L sin 6* COS0, Lsin 0sin(/), L cos6') where L = |L|, and 6 and 0 are the polar and the 
azimuthal angles of L. Note that the polar angle 6 represents the angle between the direction 
of L of the subhalo and the minor principal axis of T of its host. 

The probability density distribution of the cosines of the polar angle 6 can be obtained 
by integrating out equation (2) over L and 0 (Lee 2004): 


p(cos6*) = 




3cA. 


X 


i=l 

r*27r 


cos^ 9 


sin^ 9 cos^ 0 sin^ 9 sin^ 0 
1 -|- c — 3cAf 1 -|- c — 3cA| 1 -f c — 3 cA| 


dcj). 


( 3 ) 
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Here the polar angle 9 is forced to be in the range of [ 0 , 7 r/ 2 ] satisfying p{0) sin 9d9 = 1, 
since we care about the relative spatial orientation of the subhalo axis, but not its sign. 
Here the three Aj’s {i = 1, 2, 3) are the eigenvalues of T in a decreasing order satisfying the 
following two conditions: (i) Aj = 0; (ii) = 1. If T is a Gaussian random held 

which is true in the linear regime, one can show that Ai ~ —A 3 ~ 1/v^ and A 2 ~ 0 (Lee & 
Pen 2001). 

We adopt the following two assumptions: (i) On cluster scale, the principal axes of 
the inertia shape tensor lij of a host halo is aligned with its tidal shear tensor Tij with the 
eigenvalues being in an opposite order. In other words, the major principal axis of lij is the 
minor principal axis of Tij, and vice versa. Note that in Lee (2004), it was erroneously stated 
that it the major axis of lij is the major axis of T^j (Trujillo, Carretero, & Juncosa 2004 in 
private communication); (ii) The minor axis of a subhalo is in the direction of its angular 
momentum. 

A justihcation of the first assumption is given by the Zel’dovich approximation (Zel- 
dovich 1970) which predicts a prefect alignment between the principal axis of lij and that 
of Tij. Since the cluster-size halos are believed to be in the quasi-linear regime where the 
Zel’dovich approximation is valid, the first assumption should provide a good approximation 
to the reality. Moreover, recent N-body simulations indeed demonstrated that lij and Tij are 
quite strongly correlated (Lee & Pen 2000; Porciani, Dekel, & Hoffman 2002). Regarding the 
second assumption, there is an established theory that the spin axis of an ellipsoidal object in 
the gravitational tidal held is well correlated with its minor axis (Binney & Tremaine 1987), 
which was also conhrmed by several N-body simulations (e.g., Faltenbacher et al. 2002) 

Now that the minor principal axis of T is the major axis of the host halo, and L is 
aligned with the minor axis of the subhalo, the polar angle 6 in equation (3) actually equals 
the angle between the minor axis of the subhalo and the major axis of the host halo. Putting 
Al = A 2 = 0, and A 3 = we simplify equation (3) into 

1 I c r / 3 Al 

p{cos9) =—{1 + c)^l - - J 1c 1^1 - - sin^ 6*sin^ 0J <70, (4) 

In the asymptotic limit of c -C 1, equation (4) can be further simplihed into the following 
closed form: 

p(cos 6) = (l —sin^ 6 (5) 

Equations (4) and (5) imply that p{6) increases as 6 increases. That is, the minor axis 
of a subhalo has a strong propensity to be anti-aligned with the major axis of its host halo. 
Hence, it explains the observed alignment effect between the subhalo and the host halo major 
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axes as a consequence of the intrinsic anti-alignment between the subhalo minor axis and 
the host halo major axis. 

The value of c in equations (4) and (5) should depend on the distance from the host 
halo center (r), the subhalo mass (m) and redshift {z): c = c{r,m,z). What one can 
naturally expect is that c should decrease with r since the tidal interaction must be strongest 
in the inner part of the host halo, and that c should increase with m and r since the 
alignment effect gets reduced in the nonlinear regime. Unfortunately, it will be very difficult 
to hnd thefunctional form of c(r, ; 2 , m) as c contains all the nonlinear informations of galaxy 
evolution. We do not attempt to hnd c(r, z, m) here since it is beyond the scope of this Letter. 
Instead, we simplely assume that c is a constant, and determine the value of c empirically 
by htting equation (4) to the numerical results in §3. 


3. NUMERICAL EVIDENCES 

The data we use in this Letter is the high resolution halo simulations of Jing & Suto 
(2000). First they selected dark matter halos from their previous cosmological P^M N-body 
simulations with 256^ particles in a 100/i“^Mpc cube (Jing & Suto 1998). The halos were 
identihed using the standard friend-of-friend (FOF) algorithm, among which four halos on 
cluster-mass scales (with mass around 5 — 10 x 10 ^^/i“^Mq) were then re-simulated using 
the nested-grid P^M code which was designed to simulate high-resolution halos. The force 
resolution is typically 0.4% of the virial radius, and each halo is represented by about 2 x 
10® particles within the virial radius. We then use the SUBFIND routine of Springel et al. 
(2001) to identify the disjoint self-bound subhalos within these halos, and include those 
subhalos containing more than 10 particles in the analysis. These simulations adopted the 
“concordance” ACDM cosmology with flo = 0.3, Ua,o = 0.7, and h = 0.7. 

Using this numerical data, we hrst compute the inertia tensors as lij = Samara, 
for each host halo and its subhalos in their respective center-of-mass frames. Then, we hnd 
the directions of the eigenvectors corresponding to the largest and the smallest eigenvalues 
of the host halo and the subhalo inertia tensors, respectively, by rotating the inertia tensors 
into the principal axes frame, and determine the major axes of each host halo and the minor 
axis of its subhalos. Then, we measure the cosines of the angles, 6, between the major 
axis of the host halo and the minor axes of their subhalos by computing cos 6 = Bh ■ Sg 
where and represent the major and the minor axes of the host halo and its subhalos, 
respectively. Finally, we hnd the probability density distribution, p(cos6*), by counting the 
number density of the subhalos. When computing the probability density distribution, we 
use all the subhalos in the host halo linked by the FOF algorithm. 
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We perform the above procedure at four different redshifts: z = 0, 0.5,1 and 2. The 
total number of the subhalos Woi at each redshift is 8963, 5686, 2766, and 1469, respectively. 
Figure 1 plots the dual numerical distributions (solid dots) with the error bars. The error 
bar at each bin is nothing but the Poisson mean for the case of no alignment given as given 
as 1 /V Nbin — 1 where Win is the number of the subhalos at each bin. As one can see, 
the numerical distribution p{cosd) increases as 6 increases, revealing that the minor axes 
of substructures really tend to be anti-aligned with the major axes of their host halos, as 
predicted by the analytic model (eq.[4]) of §2. Figure 1 also plots the analytic predictions 
(solid line) and the approximation formula (dashed line) derived in §2. The horizontal dotted 
line represents the uniform distribution of cos 6 for the case of no alignment. 

We fit the analytic distributions to the numerical data points to determine the best-fit 
values of the correlation parameter c. We find c = 0.28 ± 0.01, 0.36 ± 0.02, 0.41 ± 0.02, 0.45 ± 
0.03 at z = 0, 0.5,1,1.5, respectively. The errors involved in the determination of c is given 
as the standard deviation of c for the case of no alignment effect given as Cc = a/c^/(A^ ot — 1) 
where Ntot is the number of all subhalos used to compute c. The value of c increases with 
redshifts 2 :, as expected. 

In fact the value of c = 0.3 gives quite a good fit, if not the best, not only at the present 
epoch ot z = t) but also at all earlier epochs oi z = 0.5,1,1.5, which implies that the initially 
induced anti-alignment effect is more or less conserved, reflecting the fact that the directions 
of the subhalo angular momentum are fairly well conserved. 

To understand the dependence of the alignment effect on the subhalo mass, we derive 
the same probability distribution but by using only the most massive 30 subhalos in each 
cluster {Ntot = 120 for each redshift), and find the corresponding best-£t values of c. We find 
c = 0.8 ± 0.11, 0.85 ± 0.11, 0.9 ± 0.11, 0.95 ± 0.11 at z = 0, 0.5,1,1.5, respectively. Figure 2 
plots the results. The approximation formula (eq.[5]) is excluded in this figure since the best- 
fit values of c for this case is pretty close to unity. Although the large error bars prevent us 
from making a quantitative statement, it is obvious that the anti-alignment effect is stronger 
for the case massive subhalos. 


4. SUMMARY AND DISCUSSIONS 

Although the currently popular anisotropic merging and infall scenario has provided a 
qualitative explanation for the BCG-cluster and cluster-cluster alignments (e.g., Hopkins, 
Bahcall & Bode 2005), no previous approach based on this scenario was capable of making 
a quantitative prediction for the alignment effect of cluster galaxies other than BCGs with 



the host clusters. We constructed an analytic model where the substructure alignment is a 
consequence of the tidal held of the host halo at least on the scale of cluster galaxies, and 
predicted quantitatively the strength of the alignment effect, by comparing the model with 
the results from recent high-resolution N-body simnlations. 

However, it is worth noting that we have yet to completely rule out the anisotropic 
infall model. An idealistic way to discriminate our analytic model from the anisotropic infall 
scenario would be to measure directly the correlation of the directions of the snbhalo angular 
momentum with the host halo orientations. Unfortunately, it is very difficult to determine 
the directions of the snbhalo angnlar momentnm vectors in cnrrent simulations. Because 
the rotation speed of a dark halo in simulations is only a few percent of the virial motion 
of the particles, one needs more than 10*^ particles to determine accnrately the direction 
of the angnlar momentum vector of a subhalo with an average rotation speed. In current 
simulations the subhalos have much fewer particles than 10^. This is why we used the minor 
axes of the snbhalos rather than the directions of the snbhalo angular momentum vectors to 
investigate the intrinsic alignments of snbstructures. 

Nevertheless, the strong alignment between the halo minor axes and angnlar momen¬ 
tnm vectors (Binney & Tremaine 1987) should indicate indirectly that the anti-alignments 
between the snbhalo minor axes and the host halo major axes observed in our simulations 
are likely to be caused by the host halo tidal held as our model predicts. Many N-body sim¬ 
ulations (Barnes & Efstathion 1987; Dubinski 1992; Warren et al. 1992; Bailin et ah 2005) 
have already proved that the dark matter halos rotate and their angnlar momentnm vectors 
are aligned with their minor axes. Furthermore, Bailin & Steinmetz (2004) demonstrated 
evidently that there are good internal alignments between the halo minor axes and angular 
momentum vectors measured at different radii. Therefore, the alignments between the halo 
angnlar momentum vectors and the minor axes are expected to be hold even when the outer 
parts of the halos get disrupted when they fall into larger halos as substructnres. Indeed, we 
ourselves check this effect in our simulations: we measnre the angnlar momentum vectors 
of several very massive subhalos with more than 10"^ particles, and hnd that the snbhalos 
have non-zero angular momentum and that the snbhalo minor axes are strongly aligned the 
directions of their angular momentum: the cosines of all alignment angles tnrn ont to be 
bigger than 0.6. 

Anyway, it will be dehnitely necessary to investigate directly the correlations of the 
directions of the snbhalo angular momentum with the orientations of their host halos in 
the future with larger simulation data, where the number of particles belonging to snbhalos 
should be large enough. Using larger simulation data, it will be also possible to determine the 
fnnctional form of the correlation parameter c{r,m,z). Our future work is in thisdirection. 
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Fig. 1.— Probability density distributions of the angles between the minor axes of the 
subhalo and the host halo at four different epochs ] z = 0, 0.5,1 and 1.5. In each panel, the 
dots represent the simulation results with the Poisson errors, while the solid and dashed lines 
represent the theoretical prediction (4) and the approximation formula (5), respectively. The 
horizontal dotted line corresponds to the case of no alignment effect. 
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Fig. 2.— Same as figure 1 but with the most massive 30 subhalos. 
























































